Search for: All records

Recent observations have found a growing number of hypervelocity stars with speeds of ≈1500 − 2500 km s⁻¹that could have only been produced through thermonuclear supernovae in white dwarf binaries. Most of the observed hypervelocity runaways in this class display a surprising inflated structure: their current radii are roughly an order of magnitude greater than they would have been as white dwarfs filling their Roche lobe. While many simulations exist studying the dynamical phase leading to supernova detonation in these systems, no detailed calculations of the long-term structure of the runaways have yet been performed. We used an existing AREPOhydrodynamical simulation of a supernova in a white dwarf binary as a starting point for the evolution of these stars with the one-dimensional stellar evolution code MESA. We show that the supernova shock is not energetic enough to inflate the white dwarf over timescales longer than a few thousand years, significantly shorter than the 10^5 − 6year lifetimes inferred for observed hypervelocity runaways. Although they experience a shock from a supernova less than ≈0.02 R_⊙away, our models do not experience significant interior heating, and all contract back to radii of around 0.01 R_⊙within about 10⁴years. Explaining the observed inflated states requires either an additional source of significant heating or some other physics that is not yet accounted for in the subsequent evolution. 

Abstract Type Ia supernovae (SNe Ia), critical for studying cosmic expansion, arise from thermonuclear explosions of white dwarfs, but their precise progenitor pathways remain unclear. Growing evidence supports the “double-degenerate scenario,” where two white dwarfs interact. The absence of nondegenerate companions capable of explaining the observed SN Ia rate, along with observations of hypervelocity white dwarfs, interpreted as surviving companions of such systems, provide compelling evidence for this scenario. Upcoming millihertz gravitational-wave observatories like the Laser Interferometer Space Antenna (LISA) are expected to detect thousands of double-degenerate systems, though the most compact known candidate SN Ia progenitors produce marginally detectable signals. Here, we report observations of ATLAS J1138-5139, a binary white dwarf system with an orbital period of just 28 minutes. Our analysis reveals a 1M_☉carbon–oxygen white dwarf accreting from a high-entropy helium-core white dwarf. Given its mass, the accreting carbon–oxygen white dwarf is poised to trigger a typical-luminosity SN Ia within a few million years, to evolve into a stably transferring AM Canum Venaticorum (or AM CVn) system, or undergo a merger into a massive white dwarf. ATLAS J1138-5139 provides a rare opportunity to calibrate binary evolution models by directly comparing observed orbital parameters and mass-transfer rates closer to merger than any known SN Ia progenitor. Its compact orbit ensures detectability by LISA, demonstrating the potential of millihertz gravitational-wave observatories to reveal a population of SN Ia progenitors on a Galactic scale, paving the way for multimessenger studies offering insights into the origins of these cosmologically significant explosions. 

ABSTRACT We present a detailed modelling study of CD-30°11223 (CD-30), a hot subdwarf (sdB)-white dwarf (WD) binary identified as a double detonation supernova progenitor, using the open-source stellar evolution software MESA. We focus on implementing binary evolution models carefully tuned to match the observed characteristics of the system including log g and Teff. For the first time, we account for the structure of the hydrogen envelope throughout the modelling, and find that the inclusion of element diffusion is important for matching the observed radius and temperature. We investigate the two sdB mass solutions (0.47 and 0.54 M⊙) previously proposed for this system, strongly favouring the 0.47 M⊙ solution. The WD cooling age is compared against the sdB age using our models, which suggest an sdB likely older than the WD, contrary to the standard assumption for compact sdB-WD binaries. Subsequently, we propose a possible alternate formation channel for CD-30. We also perform binary evolution modelling of the system to study various aspects such as mass transfer, orbital period evolution, and luminosity evolution. Our models confirm CD-30 as a double detonation supernova progenitor, expected to explode ≈55 Myr from now. The WD accretes an ≈0.17 M⊙ thick helium shell that causes a detonation, leaving a 0.30 M⊙ sdB ejected at ≈750 km s−1. The final 15 Myr of the system are characterized by helium accretion which dominates the system luminosity, possibly resembling an AM CVn-type system. 

Abstract Regular, automated testing is a foundational principle of modern software development. Numerous widely used continuous integration systems exist, but they are often not suitable for the unique needs of scientific simulation software. Here we describe the testing infrastructure developed for and used by the Modules for Experiments in Stellar Astrophysics (MESA) project. This system allows the computationally demanding MESA test suite to be regularly run on a heterogeneous set of computers and aggregates and displays the testing results in a form that allows for the rapid identification and diagnosis of regressions. Regularly collecting comprehensive testing data also enables longitudinal studies of the performance of the software and the properties of the models it generates. 

ABSTRACT Blue large-amplitude pulsators (BLAPs) make up a rare class of hot pulsating stars with effective temperatures of ≈30 000 K and surface gravities of 4.0–5.0 dex (cgs). The evolutionary origin and current status of BLAPs is not well understood, largely based on a lack of spectroscopic observations and no available mass constraints. However, several theoretical models have been proposed that reproduce their observed properties, including studies that identify them as pulsating helium-core pre-white dwarfs (He-core pre-WDs). We present here follow-up high-speed photometry and phase-resolved spectroscopy of one of the original 14 BLAPs, OGLE-BLAP-009, discovered during the Optical Gravitational Lensing Experiment. We aim to explore its pulsation characteristics and determine stellar properties such as mass and radius in order to test the consistency of these results with He-core pre-WD models. Using the mean atmospheric parameters found using spectroscopy, we fit a spectral energy distribution to obtain a preliminary estimate of the radius, luminosity, and mass by making use of the Gaia parallax. We then compare the consistency of these results to He-core pre-WD models generated using Modules for Experiments in Stellar Astrophysics, with predicted pulsation periods implemented using gyre. We find that our mass constraints are in agreement with a low-mass He-core pre-WD of ≈0.30 M⊙. 

Abstract Binary systems of a hot subdwarf B (sdB) star + a white dwarf (WD) with orbital periods less than 2–3 hr can come into contact due to gravitational waves and transfer mass from the sdB star to the WD before the sdB star ceases nuclear burning and contracts to become a WD. Motivated by the growing class of observed systems in this category, we study the phases of mass transfer in these systems. We find that because the residual outer hydrogen envelope accounts for a large fraction of an sdB star’s radius, sdB stars can spend a significant amount of time (∼tens of megayears) transferring this small amount of material at low rates (∼10 −10 –10 −9 M ⊙ yr −1 ) before transitioning to a phase where the bulk of their He transfers at much faster rates ( ≳10 −8 M ⊙ yr −1 ). These systems therefore spend a surprising amount of time with Roche-filling sdB donors at orbital periods longer than the range associated with He star models without an envelope. We predict that the envelope transfer phase should be detectable by searching for ellipsoidal modulation of Roche-filling objects with P orb = 30–100 minutes and T eff = 20,000–30,000 K, and that many (≥10) such systems may be found in the Galactic plane after accounting for reddening. We also argue that many of these systems may go through a phase of He transfer that matches the signatures of AM CVn systems, and that some AM CVn systems associated with young stellar populations likely descend from this channel. 

ABSTRACT We report the results from follow-up observations of two Roche-lobe filling hot subdwarf binaries with white dwarf companions predicted to have accretion discs. ZTF J213056.71+442046.5 (ZTF J2130) with a 39-min period and ZTF J205515.98+465106.5 (ZTF J2055) with a 56-min period were both discovered as subdwarf binaries with light curves that could only be explained well by including an accretion disc in their models. We performed a detailed high-resolution spectral analysis, using Keck/ESI to search for possible accretion features for both objects. We also employed polarimetric analysis using the Nordic Optical Telescope (NOT) for ZTF J2130. We did not find any signatures of an accretion disc in either object, and placed upper limits on the flux contribution and variation in degree of polarization due to the disc. Owing to the short 39-min period and availability of photometric data over 6 yr for ZTF J2130, we conducted an extensive O − C timing analysis in an attempt to look for orbital decay due to gravitational wave radiation. No such decay was detected conclusively, and a few more years of data paired with precise and consistent timing measurements were deemed necessary to constrain $$\dot{P}$$ observationally. 

Abstract Cooling white dwarfs (WDs) can yield accurate ages when theoretical cooling models fully account for the physics of the dense plasma of WD interiors. We use MESA to investigate cooling models for a set of massive and ultramassive WDs (0.9–1.3 ) for which previous models have failed to match kinematic age indicators based on Gaia DR2. We find that the WDs in this population can be explained as C/O cores experiencing unexpectedly rapid 22 Ne sedimentation in the strongly liquid interior just prior to crystallization. We propose that this rapid sedimentation is due to the formation of solid clusters of 22 Ne in the liquid C/O background plasma. We show that these heavier solid clusters sink faster than individual 22 Ne ions and enhance the sedimentation heating rate enough to dramatically slow WD cooling. MESA models including our prescription for cluster formation and sedimentation experience cooling delays of ≈4 Gyr on the WD Q branch, alleviating tension between cooling ages and kinematic ages. This same model then predicts cooling delays coinciding with crystallization of 6 Gyr or more in lower-mass WDs (0.6–0.8 ). Such delays are compatible with, and perhaps required by, observations of WD populations in the local 100 pc WD sample and the open cluster NGC 6791. These results motivate new investigations of the physics of strongly coupled C/O/Ne plasma mixtures in the strongly liquid state near crystallization and tests through comparisons with observed WD cooling.